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We investigate the dynamics of the sine-Gordon solitons perturbed by spatiotemporal external 
forces. We prove the existence of internal (shape) modes of sine-Gordon solitons when they are in 
the presence of inhomogeneous space-dependent external forces, provided some conditions (for these 
forces) hold. Additional periodic time-dependent forces can sustain oscillations of the soliton width. 
We show that, in some cases, the internal mode even can become unstable, causing the soliton to 
decay in an antisoliton and two solitons. In general, in the presence of spatiotemporal forces the 
soliton behaves as a deformable (non-rigid) object. A soliton moving in an array of inhomogeneities 
can also present sustained oscillations of its width. There are very important phenomena (like the 
soliton-antisoliton collisions) where the existence of internal modes plays a crucial role. 

PACS numbers: 05.45.-a, 05.45.Yv, 47.52.+j 
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The sine-Gordon solitons are very important in 
physics. They possess crucial applications in both parti- 
cle physics and condensed matter theory. For instance, 
in solid state physics, they describe domain walls in fer- 
romagnets, dislocations in crystals, charge density waves, 
fluxons in long Josephson junctions and Josephson trans- 
mission lines, etc. % g |, % § § 

In general, nonintegrable soliton equations (e.g. the 
equation and the double sine-Gordon) may possess 
internal degrees of freedom, which are crucial in many 
phenomena ^ ||, [| |10| . A recent discussion of inter- 
nal modes of solitary waves can be found in Ref. [p"l|| . 
However, it is well-known that the unperturbed ("pure") 
sine-Gordon equation does not have internal modes. 

A very remarkable question is the following: can ex- 
ternal forces create internal modes in the sine- Gordon 
equation? 

Recently there has been a hot debate in the scientific 
literature about the existence of internal modes of sine- 
Gordon solitons. Some authors §J, |§ H H 0> HI 
have claimed that they have found an internal quasimodc 
described as a long-lived oscillation of the width of the 
sine-Gordon soliton. 

On the other hand, a very recent and interesting paper 
is contradicting all these reports [^9) . By considering the 
response of the soliton to ac forces and initial distortions, 
Quintero et al show that neither intrinsic internal modes 
nor "quasimodes" exist in contrast to previous reports. 
We should stress that they use only time-dependent per- 
turbations in their work. 

In the present Rapid Communication we study the 



sine-Gordon equation perturbed by spatiotemporal ex- 
ternal forces: 



+ J(p t -<p xx +i 



F(x,t). 



(1) 



We will show that with some spatially inhomogeneous 
forces, the internal modes can exist for the sine-Gordon 
equation. 

We have shown in previous papers |2(], pi] ^2], |23|. |24j. 

Ba, Pq| that in equations as the following: 



(t>tt + 7<f>t ~ <t>xx + sin (j> = Fx (x) , 



(2) 



if the force Fi(x) possesses a zero x* (Fi(x*) — 0), this 
can be an equilibrium position for the soliton. If there is 
only one zero, this is a stable equilibrium position for the 

soliton if (^ dF g^ ^ > 0. For the antisoliton, it is stable 



if f^M 

V ox 



< 0. 



Let us suppose that Fi(x) is defined as: 

Fi(x) = 2(B 2 - l)smh(5a;)/ cosh 2 (fix). (3) 

This is a function with a zero in the point x* = 0. 

We have chosen this function because of the following 
properties: (i) the exact solution for the soliton resting 
on the equilibrium position can be obtained, and (ii) the 
stability problem of this soliton can be solved exactly. 
The results obtained with this function can be general- 
ized qualitatively to other systems topologically equiva- 
lent to this one. Besides, real physical systems are related 
to this example ||. For instance, in a Josephson junc- 
tion a perturbation that can be described by a function 
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of type F(x) 



dR{x) 
dx 



where R(x) is a bell-shaped func- 



tion, is an Abrikosov vortex lying in the junction's plane 
perpendicular to its local dimension [ p7| . 

A similar function can describe a local deformation of 
a charge density wave system p3| . 

Usually, function R(x) is taken as the Dirac's 6- 
function. However, if we wish to model a finite- width 
inhomogeneity, function (^) is a better choice. 

The exact stationary solution of Eq. (g), with Fi(x) 
as defined in (||), is ^ — 4arctan [exp (Bx)] . 

The stability analysis, which considers small amplitude 
oscillations around (j>k \(j>(k, x) — (f>k(x)+ f(x)e xt ] , leads 
to the eigenvalue problem [g0[ |2l], £2|, p3[ [24|, gj| p6[ : 
Lf = Tf, where L = -<9 2 + [l - 2 cosh~%Ba;)] and T = 
-A 2 - 7 A. 

This problem can be solved exactly [^9[ . The eigenval- 



ues of the discrete spectrum 
are given by the formula 




25 2t 



T„ = B 2 (A + 2An- n 2 ) - 1, 



(4) 



where A(A + 1) = 2/B 2 . 

The integer part of A, i.e. [A], yields the number of 
eigenvalues in the discrete spectrum, which correspond to 
the soliton modes (this includes the translational mode 
To, and the internal or shape modes T„ with n > |2C| , 

f^g, |§,[§, |2f, |2| . 

All this theoretical investigation produces the follow- 
ing results (note that parameter B will be our control 
parameter): For B 2 > 1, the translational mode is stable 
and there are no internal modes. If | < B 2 < 1, then 
the translational mode is unstable. However, still there 



are no internal modes 
translational mode 
ternal mode is stable 



When i < B 2 < ±, apart of the 



there is one internal mode! This in- 
In the case that B 2 < i there can 
appear many other internal modes! The exact number 

?2 For B 2 < 3 



is [A] - 1, where A(A + 1) = 2/B 



_ 5+y/l7 
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A*(A»+1)' 

the first internal mode becomes un- 



where A 
stable! 

What happens when we shift the soliton center- of-mass 
away from the equilibrium position ? 
We have the following initial problem: 

(f>(x, 0) = 4arctan{exp [B (x — xq)]} , (5) 

<Mx,0) = 0. (6) 

In the stable case (B 2 > 1) the center-of-mass of the 
soliton will make damped oscillations (for xo ^ 0) around 
the equilibrium point x = 0. 

In the case that the translational mode is unstable 
(i < B 2 < 1), the soliton will move away indefinitely 
from the equilibrium position. 

Consider the next initial problem: 

4>(x,0) = 4 arctan [exp (Bx)) + C smh(Bx) cosh~ A (Bx), 

(7) 



In this initial problem the initial soliton is deformed. 

For \ < B 2 < i we will observe oscillations of the 
soliton width (see Fig. |l|) . This is due to the fact that an 
internal mode has been excited. Eventually, due to un- 
avoidable errors in the initial conditions or to energy ex- 
change between the internal mode and the translational 
mode, the soliton will move away from the equilibrium 
position (remember that the equilibrium position is un- 
stable for the soliton center-of-mass). 




FIG. 1: Soliton's width oscillations when the internal mode 
can be excited and it is stable, i < B 2 < i 



It is important to note that the instability of the trans- 
lational mode does not mean instability of the soliton 
structure. 

We have to say that the frequency of the oscillations 
observed in the numerical simulations coincides with the 
one obtained theoretically using equation (@). The fre- 
quency of the width oscillations can be obtained using the 
equations uj\ = yTT; where Ti = i? 2 (3A— 1) — 1. All our 
experiments confirm the prediction about the frequency 
of the shape oscillations. 

The most spectacular phenomenon occurs for B 2 < 
2 a _ 5+Vrf 

A(A+1) ' 11 ~ 2 

is unstable 



&M) = o. 



In this case, the first internal mode 
If we study the evolution of the soliton from 
the initial conditions (0)-(^) we will observe the destruc- 
tion of the soliton (see Fig. |^). Two solitons move away 
(in different directions) to "infinity" (or to the bound- 
aries of the system and an antisoliton is formed in the 
place of the original soliton remaining there stabilized. 
In fact, the condition B 2 < 1 implies stability for the 
center-of-mass of an antisoliton. 

Note that in these situations, the sine-Gordon solitons 
do not behave as rigid objects, which is what is expected 
from them in general ^TJ . The initial distortions of the 
width of the soliton will eventually be damped due to 
dissipation. 

Once the internal modes are possible, as in Eqs. (||)- 
(|) with i < B 2 < |, we need time-dependent external 
forces to sustain the oscillations of the soliton width. 

On the other hand, if we wish the soliton to remain in 
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FIG. 2: Soliton's destruction when the internal mode is un- 



stable, B 2 < 2 , where A, = 8±2j3I. 



FIG. 3: Soliton profiles for different time instants governed 
by Eq. (§) (S = 0.5, D = 0.2, 7 = 0.1, f = 0.7, u = 
0.55, E = 0.7). 



some spatially localized zone, we need stable equilibrium 
positions for the center-of-mass of the soliton. 

Let us consider the following spatiotemporal perturba- 
tion: 



(f>tt + 7& " 4>*x + sin0 = F 2 (x) + F 3 (x, t), (9) 



where 



F 2 (x) 



Fi(x) 



A 



cosh [B (x + xi)] 
D A 



cosh [B (x — xi)] 



, if — X\ < x < x\, 
D , if x < —x±, 

, if x > Xi, 



and F 3 (x,t) = f cos(wi)[l/cosh 2 (_E(x + xi)) + 
l/cosh 2 (S( x — Xi)). The space-dependent force F 2 (x) 
creates a double-well potential for the soliton. At the 
same time, in the interval — x\ < x < Xi, we have the 
same force F%(x), which was sufficient for the existence 
of the internal mode. 

Actually, other forces can be used to excite the internal 
mode. In fact, if we have a function F(x) that can mimic 
approximately the behavior of function F% (x) (specially 
in the interval — xi < x < xi, where — x\ and x\ are the 
extrema of function -Fi(x)) when B satisfies the condition 
g < B 2 < ~, then this function is good for exciting the 
internal mode. And note that the behavior of function 
-F\(x) in the interval —x% < x < x\ is a very common 
behavior for a function in an interval where there is a 
zero and two extrema. 

The time-dependent force F^(x,t) will cause the soli- 
ton width to oscillate. The center-of-mass of the soliton 
will also oscillate, jumping between the potential wells 
created by force F 2 (x). 

Although the soliton is not always in the interval 
— x\ < x < Xi, it will return to this interval regularly. 
While the soliton is in this interval, all the conditions 
hold for the internal mode to be excited. 



Figure]^ shows the extraordinary deformations suffered 
by the soliton in some cases. 

However, we have corroborated that other spatiotem- 
poral forces also can sustain the soliton width oscilla- 
tions. This includes the "ubiquitous" force F-^it) — 
/ocos(wi). 

The propagation of solitons in disordered media has 
been studied intensively in recent years (|(], Q . 
Consider the equation 



4>tt + l4>t - 4>xx + sin0 = F(x). 



(10) 



where F(x) is defined in such a way that it possesses 
many zeroes, maxima and minima. This system describes 
an array of inhomogeneities. 

For our study, we have defined F(x) in the following 
way: 



F(x) 



^ 4 (1 - -B 2 ) — (11) 



n— — q 



^ e 2B(x+x n ) + iy 



where x n = (n + 2) log (V2 + l) /B (n = —q,—q + 
1 • • • , q— 1, q), and q + 2 is the number of extrema points 
ofF(x). 

In our array, there is a "superposition" of the "disor- 
der" with a dc component which will cause the soliton to 
move to the right all the time. When the soliton is mov- 
ing over intervals where < 0, the internal mode can 
be excited. In fact, the points Xi where F(xi) = and 
dF ^^ < 0, are "barriers" which the soliton can over- 
come due to its kinetic energy. These "collisions" with 
the barriers will excite the internal modes if in these in- 
tervals the function F(x) mimics the behavior of -Fi(x) 
when i < B 2 < i The simulations show that the width 
of the soliton will perform sustained oscillations during 
its motion in a disordered medium. 

We have shown that in sine-Gordon equations per- 
turbed by inhomogeneous (space-dependent) forces F(x), 
the solitons can possess internal modes. 
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Some of our results are in agreement with previous 
works (H, [p. In fact, in Eq. (§), with Fi(a;) as in §), 
if we put B 2 = 1, then there are no external forces, and 
fromEq. (|) we obtain that there are no internal modes 
in that case neither. 

Moreover, even when there is an inhomogcneous exter- 
nal force, not for every F(x) we have internal modes. 

For instance, if F(x) has a zero which corresponds to 
a stable equilibrium position for the soliton, even then 
the internal modes are impossible. This explains why it 
has been so difficult to find sine-Gordon internal modes. 
For the existence of internal modes for sine-Gordon soli- 
tons we need zeroes of function F(x) that corresponds to 
unstable positions for the center-of-mass of the soliton. 
When the soliton center-of-mass is very close to an unsta- 
ble equilibrium position, there is a pair of forces acting 
in opposite directions on the "body" of the soliton. This 
pair of forces should be sufficiently large to stretch the 
soliton "body" , such that the soliton internal mode can 
be excited. 

Function F{x) can possess several zeroes correspond- 
ing to unstable and stable equilibrium positions. For 
instance, we have studied a force F(x) that creates a 
double-well potential for the soliton. 

Periodic time-dependent forces (besides the spatially 
inhomogeneous forces) can sustain the oscillations of the 



soliton width. 

A soliton moving in an array of inhomogeneities can 
also undergo sustained oscillations of its width. 

All this is possible because the internal mode of the 
soliton can exist when it is moving in media where there 
are inhomogeneous space-dependent forces with unstable 
equilibrium positions. 

Nonetheless, we have discovered another more remark- 
able phenomenon: The sine-Gordon internal mode not 
only can exist for some external forces, but (in some sit- 
uations) it can become unstable. If we have an unstable 
equilibrium position for the soliton center-of-mass and 
the pair of forces acting on the soliton is too large, then 
the soliton can be destroyed. 

When the soliton is destroyed, it can be transformed 
into an antisoliton and two new solitons. The topologi- 
cal charge is conserved. We had found this phenomenon 
before for the </> 4 -equation pll |22j. However, here for the 
first time we have shown not only that the sine-Gordon 
soliton internal mode can exist, but that it can become 
unstable and destroy the soliton. This is a spectacular 
manifestation of the fact that the sine-Gordon soliton can 
behave as a deformablc (non-rigid) object. 

A. Bellorm would like to thank CDCH-UCV for sup- 
port under project PI-03-11-4647-2000. 
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